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Abstract 

In this paper, based on the works of Capozziello et al., we have studied the Noether symmetry 
approach in the cosmological model with scalar and gauge fields proposed recently by Soda et al. 
The correct Noether symmetries and related Lie algebra are given according to the minisuperspace 
quantum cosmological model. The Wheeler-De Witt (WDW) equation is presented after quanti¬ 
zation and the classical trajectories are then obtained in the semi-classical limit. The oscillating 
features of the wave function in the cosmic evolution recover the so-called Hartle criterion, and the 
selection rule in minisuperspace quantum cosmology is strengthened. Then we have realized now 
the proposition that Noether symmetries select classical universes. 

PACS numbers: 98.80.Cq, 98.80.Qc, 04.50.Kd, 95.36,+x 


1 


The Noether symmetry approach 


Q 


in general re’ 

n 

ment of the crucial role of symmetries in physics [2- 


ativity and cosmology is the endorse- 
6]. As pointed out by Capozzicllo et 


al. ^-4], such an approach can provide fundamental consideration on the related models 
in terms of symmetries, and a selection rule J 4 ] will come into being to recover the classi¬ 
cal behaviors of cosmic evolution due to the oscillating features of the cosmological wave 
function. That is, the correlated region in the configuration space of dynamical variables 
is selected and a classical observable universe emerges. In the minisuperspace description, 
the so-called Hartle criterion [7] is applied to the solutions of the Wheeler-De Witt (WDW) 
equation and selects the classical trajectories Q. 

Moreover, the vector field (or gauge field) has been taken into account in cosmology, es¬ 


pecially for the early universe [8 


10]. In the cosmological models with the vector field(s), the 


anisotropic hair, perturbation, non-Gaussianity and other aspects are studied extensively. 
Furthermore, if we are in the case of gauge field, that is, the gauge symmetry is conserved 
in the system, as proposed by Soda et al. 8-Ty], the exact anisotropic power-law solutions 
have been obtained when both the potential function for the scalar field and the gauge ki¬ 
netic function are exponential type, and an attractor arises for a large parameter region [ 9 ]. 
The anisotropic hair related to the so-called cosmic no-hair theorem and the mechanism of 
magnetogenesis can also be discussed on this stage. 

Furthermore, the Lie point symmetries of second-order partial differential equations, in¬ 
cluding the Schrodinger and the Klein-Gordon equations have been studied in a general 
Riemannian space by Paliathanasis et al. [ll, 12], and these symmetries are related to the 
Noether point symmetries of the classical Lagrangian. By employing the Lie symmetries 
of the WDW equation, a general family of hyperbolic scalar field potentials is also dis¬ 
cussed within the framework of perfect fluid in Ref. [ll], and the exact solutions of the field 
equations could be obtained based on the Lie symmetries of the WDW equation. 

Then, if we want to discuss the Noether approach in cosmological model with scalar and 

9] 


vector fields, we should start with the following Lagrangian 

1 _ . 1 


S — J d A xyf^g 


M 2 

-fR - -(dMW) - K(0) - 


( 1 ) 


where the flat Robertson-Walker background geometry is assumed, g, R and a(t) are the 
determinant of the metric, the Ricci scalar and the scale factor respectively, the reduced 
Planck mass takes M p = 1, while V(4>) is a potential of the scalar field </> and /(</>) is 
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the gauge kinetic function, a coupling between the scalar held and the strength tensor 
F /tl/ = d fJ A v — d v A p for the vector held A p . Such a model has been studied recently in 
Noether symmetry approach by Vakili in Ref. js| - Here, based on the works of Capozziello 
et al., we would like to begin anew to provide the correct Noether symmetries and related Lie 
algebra thanks to the minisuperspace quantum cosmological model. Then we try to obtain 
the solutions of the Wheeler-De Witt (WDW) equation and present the classical trajectories 
in the semi-classical limit. The Hartle criterion should be applied and the selection rule is 
reinforced. 

As in Refs, jg, 9], we take the vector held A p along the x-axis A p = (0, A x (t), 0, 0), and 
the scalar held 0 = 0(t), both of them are considered as homogeneous helds. After some 
manipulations and do an integration by parts, a point-like Lagrangian is obtained in the 
following form 

£ = —3 aa 2 + -a 3 0 2 + - a/ 2 (0)A 2 — a 3 H(0). (2) 

The configuration space for such a Lagrangian is Q = (a, 0, A), and then cosmological 
dynamics can be reached on such a three dimensional FRW minisuperspace. 

According to Capozziello et al. [2|, the vector held X 


X = cd(g)A + cP(g)A 
y ' dq 1 y ' dq l 


( 3 ) 


could be acted on a Lagrangian £ dehned on the tangent space of configurations TQ, then 
the Lie derivative of £ is 


L X C = XC = a‘fe)0 + a‘(q)AL. 

If LxjC = 0, we call X a symmetry for the dynamics of £. ft is clear that £ in Eq. 
not depend on A, which is then cyclic. So we have the symmetry definitely 

X = — 

3 dA' 

that is, the related a 3 is a constant and independent of variables a, 0 and A. 
Considering the model Eq. (J2J) presented here, we have 


( 4 ) 

does 

( 5 ) 


id 


d 


d 


m3 


d 


d 


X — (X — -b (X ~x ~~,—b (X 7———b (X 7——b (X —;—b OL - 


( 6 ) 


da dcp dA da d<p dA 
If the Lie derivative of £ vanishes along X, that is, LxjC = 0 is satished, a constant of 
motion can be deduced and the Noether theorem holds. 
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After a straightforward calculation, we get a quadratic expression of a, 0 and A. Then 
each coefficient should be zero and a system of partial differential equations for a 1 {i = 1,2,3) 


is obtained. 


a 1 + 2a-—— = 0, 
da 

3a 1 + 2a—— = 0, 
dtp 


da 1 


,da 2 


dtp da 
da 1 „ 9/ <9a 3 

- 6 m +/w &t = 0 ' 

9 <9a 2 „ 9/ , N <9a 3 

“ M + /w 7m= 0 ’ 


a 1 /^) + 2 aa 2 f\(f)) + 2a/(0) 


<9a 3 

a4 


= 0 , 


3aV(0) + aaV(0) = 0, 


(7) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 


where the prime is a derivative with respect to the scalar field 0. 

To recover the symmetry Eq. (J5|). it requires that a 3 is independent of variables a, 0 and 
A. Then from Eq. (flOl) . da 1 jdA = 0, and due to Eq. (0), 


a 1 (a, 0) = a ^ 2 F(0), 


(14) 


where F is an arbitrary function dependent of 0 only. From Eq. ([IT]) , do? jdA 
to Eq. (JHJ), 

ct 2 (a, 0) = — ^a -3 / 2 J F(0)d0. 

From Eq. (TT2]h a 1 / + 2aa 2 f = 0, and 

f a 1 

f 2 aa 2 

From Eq. (1X3]) . 3adX + aa 2 V' = 0, and 

V' _ 3a 1 

V ~ ~aa 7 ' 


0, and due 


(15) 


(16) 


(17) 


Then from Eq. (f9[), 



we can obtain the following solution for F(0) 


0, 


(18) 


F(4>) = b ie “* + b 2 e~^, 


(19) 
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where b\ and b 2 are constants of integration, and a; 2 = 3/8. The expressions for a 1 and a 2 
are also found 

(he"* + b 2 e~^) , (20) 


a 1 = a -1/2 (h,*** 


or 


= a~ 3/2 (6ie“ 0 - & 2 e“ w *) . 

ZUJ 


( 21 ) 


It should be emphasized that there is no A in the expressions of a 1 and a 2 , and the symmetry 
Eq. (j5J) is recovered. From Eq. (fT51) 

/' 1 b 1 e^ + b 2 e-^ 


then the solution of /(</>) is 


From Eq. (U71) 


then the solution of V(<j>) is 


/ 8a; bie w< ^ — b 2 e~ u ^ ’ 

/(0) = (fe ie ^-6 2 e-^) 1/3 

V' _ 3 b ie ^ + b 2 e~ 

V 4o; b\e UJ( t > — & 2 e - ^ ’ 


V((/>) = (h - b 2 e-“*) 2 . 
So we reach a basis of symmetries on T Q, 

X , = a 'X K + X a -m e -^JL + (,, 


X 2 = - X a -W e ^l_ + (, 

aa 2a; ^ 


a 1 e 


X = — 
3 &4‘ 


The corresponding constants of motion associated with the above symmetries are 

Ki = L'V “*a - ua- 3 ^e-“U, 

K 2 = 5 a i/2 e "* 4 + W a 3/2 e"*^, 

A' 3 = a/ 2 A 


( 22 ) 


(23) 

(24) 

(25) 


XrX 1 

^A a -3/2 e V) 

^2a; ) 

■ a 

(26) 

1 

^ 1^ 

3 

(J-a~ 3/2 e «A 

• a 
a? 

(27) 




(28) 


(29) 

(30) 

(31) 


It can be shown that all the symmetries are commuted with each other, that is, the 
following Lie algebra is satisfied, 


[Xi, Xj] = 0, i,j = 1,2,3. 


(32) 
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And the corresponding constants of motion also close the same algebra in terms of Poisson 
bracket, 

[Ki,Kj]= 0, i,j = 1,2,3. (33) 

According to Capozzicllo et ah jsj, to identify the cyclic variables, we make a transfor¬ 
mation on the vector held Eq. (J6J) (or Eq. f[26jl ) for X\, 


w = a 3/2 e^, 


z = a 3/ V^, 


A = A, 


(34) 


to satished that 

ix 1 dw = 3, ix 1 dz = 0, ix^A = 0. (35) 

The configuration space is transformed from Q = (a, </>, A) to Q = ( w , z, A) correspondingly. 
Then the Lagrangian Eq. ([2]) transformed into 


C = -\wz ~Z 2 + h 2/3 i 2 , 

o Z 

and the conjugated momenta of dynamic variables w, z and A are respectively 


7T w 


dC 

dw 




3 ’ 


dC 

dz 

dC 


= dl = ~3 W > 


7T A = ^ = Z^A. 

8A 

In the Hamiltonian dynamics, we should get the Hamiltonian firstly 

U = -^7 T z + ^Z~ 2/3 7T 2 a + Z 2 , 

and due to the Noether symmetries, the constants of motion are 


77 w ^0 j 


7 Ta ~ Si- 


We can quantize the system as follows, 


7T W ^ 7 Tuu 


7T 2 -> Ttz = -id z , 


77,4 -t Tt A = ~idA, 


H —>■ 'H(q j , —id q j), 

then the Wheeler-De Witt (WDW) equation could be obtained 

+ ^z~ 2/3 (-id A ) 2 + z 2 | *) = 0. 


(36) 

(37) 

(38) 

(39) 

(40) 

(41) 

(42) 

(43) 

(44) 
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Separating the variables via the constants of motion, 


*> =1 V(w)) I x(A)) I Z(z)) cx | £(z)), (45) 


the solution of £(z) is found 


£(z) = exp 


2Ef 


,1/3 


+ i 


9Sr 


(46) 


Such a wave function with the oscillating feature recovers the so-called Hartle criterion [7] 
and allows us to get a classical observable universe. 

If we now take the action (or the Hamiltonian principal function) as 

S = E 0 u + ^z V 3 + -!-z\ (47) 

the corresponding momenta which conjugate the generalized coordinates are obtained 


_ dS_ 
dw 


£ 


o, 


K z 


dS_ 

dz 


2Ef 

3 £^ 

dS 



2/3 + 



(48) 

(49) 

(50) 


Separating the variables via above formulae, we can solve Hamilton-Jacobi (HJ) equation in 
the semi-classical limit. The classical trajectories in the configuration space, Q = (w,z,A), 
then can be expressed 


z{t) = fcit + k 2 , 

(51) 

A(t ) = c{k\t + k 2 ) 1/3 + H 0 , 

(52) 

Ci + c 2 (k\t + k 2 ) l d 3 + c 3 (/cit + k 2 ) 3 , 

(53) 


where k 2 , A 0 and C\ are constants of integration, while ki, c, c 2 and c 3 are dependent of 
constants of motion £ 0 and Ex. If we go back to Q = (a,4>,A), we can get the classical 
cosmological solutions, 


a (t) — [di(t — t 0 ) + d 2 (t — t 0 ) 4//3 + d 3 (t — t 0 ) 4 ] ^ , (54) 

m = --In-ht+h- (55) 

u [di(i - to) + d 2 (t - t 0 ) 4/3 + d 3 (t - to) 4 ] 7 

where to and di (i = 1,2,3) could be looked on as constants of integration and should 
be determined by observations. Then the oscillating regime is selected and the classical 
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behaviors are recovered. One of the main results is that the late time acceleration could 
be obtained from the classical solution, and the roll of scalar held may be more important 
than that of vector held as discussed in Ref. [5]. As to the Theorem in Ref. [4], Noether 
symmetries select classical universes. Such a picture has now been realized here in the 
cosmological model with vector held(s) proposed by Soda et al AS- 

Moreover, the approach given here can be generalized naturally and consistently to the 


case of phantom held 


13] and other cases (such as the so-called quintom held) rather than 


other approaches. If we take the following Lagrangian, 


S — J d A x\J—g 


M 2 1 i 

-fR - e-(d^Xd'V) - R(0) - 


(56) 


then e = +1 is the standard case just discussed here, while e = — 1 is the phantom held case 
ji0. Same approach could be applied naturally though different behaviors are expected, 
and related energy conditions are also to be discussed. Furthermore, the related WDW 
equation is more complex, then we shall discuss the Lie symmetries of the WDW equation 
as those in Ref. [hi] to identify the properties of the solution. All these should be prepared 
systematically together with various other aspects and will be placed somewhere else. 

In summary, we have realized the proposition that Noether symmetries select classical 
universes in the cosmological model with vector (gauge) held(s). Based on the correct 
Noether symmetries, we give the related Lie algebras for symmetries and constants of motion 
according to the minisuperspace quantum cosmological model. We solve the Wheeler-De 
Witt (WDW) equation after quantization and select a subset of the solutions with oscillating 
behaviors. Then the classical observable universe is recovered. 


Finally, I would like to give a brief comment on arXiv: 1410.3131 [5]. 

Though it is an interesting work, frankly, the symmetries in arXiv:1410.3131 are incorrect. 

Firstly and evidently, -A is not a symmetry there, but in fact it is. While the constant 
of motion related to this symmetry, called P$a or Pa, is used in the calculation and plays a 
crucial role to obtain the final results and conclusions. It is contradictive. We know that a 
constant of motion should be corresponded to a symmetry. If there is no such a symmetry, 
we cannot get the related constant of motion logically. 

Secondly, other symmetries are also not correct, and no related constants of motion can 
be given. The follow-up calculations seem not to be directly based on these constants of 
motion according to the method used in the paper, the incorrectness looks not so evident 








like the one above. So the Lie Algebra is not given, in fact cannot be given there, nor the 
WDW equation and so on. 

Importantly, if the constant of motion used in the paper, called P 0A or P A , could be 
obtained, it should be satisfied that, based on its own derivation, 

7 = la T oa-y^/ 3 T ah*"*/ 3 = 0, (57) 

that is, 

— “ ± u(j) = 0. (58) 

z a 

While the one in the paper, eq. (38), is (the quantity Q is set to be zero) 

- ±uJ<j) + 3^ = 0. (59) 

2 a A 

It is clearly that ^ ^ 0, so no any result like P$ A or P A given in the paper can be obtained. 

The above equation is crucial to obtain the final results. Without such an equation, no 
any of the final results shown in the paper could be obtained. Taking some specific values for 
the constants, the equations of the system become solvable and the results look like correct 
ones. The results just happen to be superficially, the calculations are incorrect. 

The correct symmetries and the related Lie algebra are given in the present paper, the 
cosmological variables are calculated by a different method based on the works of Capozziello 
et al. 

Acknowledgements 

The author is grateful to Professor Fan Wang, Dr. Wei Chen, Dr. Yu Jiang, Dr. Feng- 
Yao Hou, Dr. Lei Feng, Dr. Yu-Peng Yang, Dr. Yan Yan, Dr. Peng Dong for useful 
discussion and also to Mr. Yi-Qiao Dong and Mr. K.B. Kim. This work is supported in 
part by the National Natural Science Foundation of China (Grant Nos. 11373068, 10973039 
and 10447114). 


[1] V.I. Arnold, Mathematical Methods of Classical Mechanics, translated by K. Vogtmann and 
A. Weinstein, 2nd Edition, Springer-Verlag, New York, 1989. 


9 



[2] S. Capozziello, M. De Laurentis, arXiv: 1308.1208 [gr-qc]; 

S. Capozziello, A.N. Makarenko, S.D. Odintsov, Phys. Rev. D 87 (2013) 084037, 
arXiv: 1302.0093 [gr-qc]; 

S. Capozziello, M. De Laurentis, S.D. Odintsov, Eur. Phys. J. C 72 (2012) 2068, 
arXiv: 1206.4842 [gr-qc]; 

S. Capozziello, M. De Laurentis, Phys. Rept. 509 (2011) 167, arXiv:1108.6266 [gr-qc]; 

S. Capozziello, A. Stabile, A. Troisi, Class. Quantum Gravity 24 (2007) 2153, 
arXiv:gr-qc/0703067; 

S. Capozziello, G. Marmo, C. Rubano, P. Scudellaro, Int. J. Mod. Phys. D 6 (1997) 491, 
arXiv:gr-qc/9606050. 

[3] A. Borowiec, S. Capozziello, M. De Laurentis, F.S.N. Lobo, A. Paliathanasis, M. Paolella, A. 
Wojnar, arXiv: 1407.4313 [gr-qc]; 

S. Capozziello, M. De Laurentis, S.D. Odintsov, arXiv: 1406.5652 [gr-qc], to appear in Mod. 
Phys. Lett. A; 

A. Paliathanasis, M. Tsamparlis, S. Basilakos, S. Capozziello, Phys. Rev. D 89 (2014) 063532, 
arXiv: 1403.0332 [astro-ph.CO]; 

A. Paliathanasis, S. Basilakos, E.N. Saridakis, S. Capozziello, K. Atazadeh, F. Darabi, M. 
Tsamparlis, Phys. Rev. D 89 (2014) 104042, arXiv:1402.5935 [gr-qc]; 

S. Basilakos, S. Capozziello, M. De Laurentis, A. Paliathanasis, M. Tsamparlis, Phys. Rev. D 
, arXiv:1311.2173 [gr-qc]; 

S. Capozziello, N. Frusciante, D. Vernieri, Gen. Relativ. Gravit. 44 (2012) 1881, 
arXiv: 1204.4650 [gr-qc]; 

S. Capozziello, P. Martin-Moruno, C. Rubano, AIP Conf. Proc. 1122 (2009) 213, 
arXiv:0812.2138 [gr-qc]; 

H. Motavali, S. Capozziello, M. Rowshan Almeh Jog, Phys. Lett. B 666 (2008) 10, 
arXiv:0807.0347 [gr-qc]; 

S. Capozziello, A. De Felice, J. Cosmol. Astropart. Phys. 0808 (2008) 016, arXiv:0804.2163 
[gr-qc]; 

S. Capozziello, V.I. Man’ko, G. Marmo, C. Stornaiolo, Gen. Relativ. Gravit. 40 (2008) 2627- 
2647, arXiv:0706.3018 [gr-qc]; 

S. Capozziello, S. Nesseris, L. Perivolaropoulos, J. Cosmol. Astropart. Phys. 0712 (2007) 009, 


10 




arXiv:0705.3586 [astro-ph]; 

S. Capozziello, G. Lambiase, Gen. Relativ. Gravit. 32 (2000) 295, arXiv:gr-qc/9912084; 

S. Capozziello, R. de Ritis, C. Rubano, P. Scudellaro, La Rivista del Nuovo Cimento 19/4 
(1996) 1; 

S. Capozziello, R. de Ritis, Class. Quantum Grav. 11 (1994) 107; 

S. Capozziello, R. de Ritis, Nuovo Cim. B109 (1994) 795; 

S. Capozziello, R. de Ritis, Nuovo Cim. B109 (1994) 783; 

S. Capozziello, R. de Ritis, Phys. Lett. A 177 (1993) 1; 

S. Capozziello, R. de Ritis, Int. J. Mod. Phys. D 2 (1993) 373; 

M. Demianski, R. De Ritis, C. Rubano, P. Scudellaro, C. S tornaiolo, Phys. Rev. D 46 (1992) 
1391; 

R. de Ritis, G. Marmo, G. Platania, C. Rubano, P. Scudellaro, C. Stornaiolo, Phys. Rev. D 
42 (1990) 1091; 

R. de Ritis, G. Marmo, G. Platania, C. Rubano, P. Scudellaro, C. Stornaiolo, Phys. Lett. A 
149 (1990) 79. 

[4] S. Capozziello, G. Lambiase, Gen. Relativ. Gravit. 32 (2000) 673, arXiv:gr-qc/9912083, 

[5] B. Vakili, Phys. Lett. B 738 (2014) 488. 

[6] T. Christodoulakis, N. Dimakis, Petros A. Terzis, B. Vakili, E. Melas, Th. Gram-menos, Phys. 
Rev. D 89 (2014) 044031, arXiv:1309.6106 [gr-qc]; 

T. Christodoulakis, N. Dimakis, Petros A. Terzis, G. Doulis, Th. Grammenos, E. Melas, A. 
Spanou, J. Geom. Phys. 71 (2013) 127, arXiv: 1208.0462 [gr-qc]; 

B. Vakili, F. Khazaie, Class. Quantum Gravity 29 (2012) 035015, arXiv: 1109.3352 [gr-qc]; 

B. Vakili, Int. J. Theor. Phys. 51 (2012) 133, arXiv:1102.1682 [gr-qc]; 

B. Vakili, Phys. Lett. B 664 (2008) 16, arXiv:0804.3449 [gr-qc]; 

F. Darabi, K. Atazadeh, A. Rezaei-Aghdam, Eur. Phys. J. C 73 (2013) 2657, arXiv: 1304.2926 
[qr-qc]. 

[7] J.B. Hartle, in Gravitation in Astrophysics, Gargese 1986, eds. B. Carter, J.B. Hartle, Plenum, 
New York, 1986; 

J.B. Hartle, in Gravitation and Quantization, Proceeding of the 1992 Les Houches Summer 
School, Eds. B. Julia, J. Zinn-Justin, Les Houches Summer School Proceedings, vol. LVII, 
North-Holland, Amsterdam, 1995, arXiv:gr-qc/9304006; 


11 



D. Craig, J.B. Hartle, Phys. Rev. D69 (2004) 123525, arXiv:gr-qc/0309117; 

J.J. Halliwell, J.B. Hartle, Phys. Rev. D 41 (1990) 1815; 

J.J. Halliwell, Nucl. Phys. B 266 (1986) 228; 

J.J. Halliwell, Nucl. Phys. B 286 (1987) 729; 

J.J. Halliwell, Phys. Rev. D 36 (1987) 3626; 

J. J. Halliwell, in Quantum Cosmology and Baby Universes, Eds. S. Coleman, J.B. Hartle, T. 
Piran, S. Weinberg, World Scientific, Singapore, 1991; 

T.P. Singh, T. Padmanabhan, Ann. Phys. (N.Y.) 196 (1989) 296. 

[8] S. Kanno, J. Soda, M-a. Watanabe, JCAP 0912 (2009) 009, arXiv:0908.3509 [astro-ph.CO]; 
M-a. Watanabe, S. Kanno, J. Soda, Phys. Rev. Lett. 102 (2009) 191302, arXiv:0902.2833 
[hep-th]; 

A. Maleknejad, M.M. Sheikh-Jabbari, J. Soda, Phys. Rept. 528 (2013) 161, arXiv:1212.2921 
[hep-th]. 

[9] S. Kanno, J. Soda, M-a. Watanabe, JCAP 1012 (2010) 024, arXiv: 1010.5307 [hep-th], 

[10] M-a. Watanabe, S. Kanno, J. Soda, Prog. Theor. Phys. 123 (2010) 1041, arXiv: 1003.0056 
[astro-ph.CO]; 

K. Yamamoto, M-a. Watanabe, J. Soda, Class. Quantum Grav. 29 (2012) 145008, 
arXiv: 1201.5309 [hep-th]; 

K. Murata, J. Soda, JCAP 1106 (2011) 037, arXiv: 1103.6164 [hep-th], 

[11] A. Paliathanasis, M. Tsamparlis, S. Basilakos, J.D. Barrow, arXiv: 1503.05750 [gr-qc]. 

[12] M. Tsamparlis, A. Paliathanasis, J. Phys. Conf. Ser.: Proceedings of GADEIS VII, Larnaca, 
Cyprus, June 2014, arXiv:1501.06047 [math.AP]; 

A. Paliathanasis, Symmetries of Differential equations and Applications in Relativistic Physics, 
PhD Thesis, University of Athens (2014), arXiv:1501.05129 [gr-qc]; 

A. Paliathanasis, M. Tsamparlis, M.T. Mustafa, Int. J. Geom. Methods Mod. Phys. 12 (2015) 
1550033, arXiv:1411.0398 [math-ph]; 

M. Tsamparlis, A. Paliathanasis, Ashgar Qadir, Int. J. Geom. Methods Mod. Phys. 12 (2015) 
1550003, arXiv: 1407.4601 [math.AP]; 

A. Paliathanasis, M. Tsamparlis, Int. J. Geom. Methods Mod. Phys. 11 (2014) 1450037, 
arXiv: 1312.3942 [math-ph], 

[13] S. Capozziello, E. Piedipalumbo, C. Rubano, P. Scudellaro, Phys. Rev. D 80 (2009) 104030, 


12 


arXiv:0908.2362 [astro-ph.CO]. 

[14] R.R. Caldwell, Phys. Lett. B 545 (2002) 23, arXiv:astro-ph/9908168. 


13 



